On preconditioned eigensolvers and Invert-Lanczos processes
نویسنده
چکیده
This paper deals with the convergence analysis of various preconditioned iterations to compute the smallest eigenvalue of a discretized self-adjoint and elliptic partial differential operator. For these eigenproblems several preconditioned iterative solvers are known, but unfortunately, the convergence theory for some of these solvers is not very well understood. The aim is to show that preconditioned eigensolvers (like the preconditioned steepest descent iteration (PSD) and the locally optimal preconditioned conjugate gradient method (LOPCG)) can be interpreted as truncated approximate Krylov subspace iterations. In the limit of preconditioning with the exact inverse of the system matrix (such preconditioning can be approximated by multiple steps of a preconditioned linear solver) the iterations behave like Invert-Lanczos processes for which convergence estimates are derived.
منابع مشابه
A comparison of algorithms for modal analysis in the absence of a sparse direct method
The goal of our report is to compare a number of algorithms for computing a large number of eigenvectors of the generalized eigenvalue problem arising from a modal analysis of elastic structures using preconditioned iterative methods. The shift-invert Lanczos algorithm has emerged as the workhorse for the solution of this generalized eigenvalue problem; however a sparse direct method is require...
متن کاملState-of-the-art eigensolvers for electronic structure calculations of large scale nano-systems
The band edge states determine optical and electronic properties of semiconductor nano-structures which can be computed from an interior eigenproblem. We study the reliability and performance of state-of-the-art iterative eigensolvers on large quantum dots and wires, focusing on variants of preconditioned CG, Lanczos, and Davidson methods. One Davidson variant, the GD + k (Olsen) method, is ide...
متن کاملA Comparison of Numerical Implementations of the Eigenstate Expansion Method for Quantum Molecular Dynamics Simulations
We investigate the efficient computation of a few of the lowest eigenvalues of a symmetric eigenvalue problem occurring in quantum dynamical molecular simulations. The large sparse system of order n is highly structured such that its multiplication with a vector costs O(n logn) floating point operations only. We compare a number of eigensolvers: subspace iteration, two variants of the restarted...
متن کاملModern Preconditioned Eigensolvers for Spectral Image Segmentation and Graph Bisection Workshop on Clustering Large Data Sets Third IEEE International Conference on Data Mining (ICDM 2003)
Known spectral methods for graph bipartition and image segmentation require numerical solution of eigenvalue problems with the graph Laplacian. We discuss several modern preconditioned eigenvalue solvers for computing the Fiedler vectors of large scale eigenvalue problems. The ultimate goal is to find a method with a linear complexity, i.e. a method with computational costs that scale linearly ...
متن کاملFast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion
In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the...
متن کامل